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Interdependencies between experimental spectra, representing line or plane projections of elec- 
tronic densities, are derived from their consistency and symmetry conditions. Some additional 
relations for plane projections are obtained by treating them as line projections of line projections. 

The knowledge of these dependencies can be utilised both for an improvement of experimental 
data and for a verification of various techniques used for correcting e.g. two-dimensional (or one- 
dimensional) angular correlation of annihilation radiation spectra and Compton scattering profiles. 
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I. INTRODUCTION 

All spectra, representing projections of the same den- 
sity, must be interdependent. This so-called consistency 
condition (CC) has been considered for reconstructing 
densities p(p) from their line projections [1-3]. In the case 
of angular correlation of annihilation radiation (ACAR) 
spectra or Compton profiles ( CP) measurements this con- 
dition is satisfied if spectra are measured up to such a 
momentum p max above which p(p) is isotropic and all 
projections have the same values. 

The CC is automatically imposed on the experimen- 
tal data via the reconstruction of p(p). However, it can 
be utilised for checking (before the reconstruction) if data 
were measured and next corrected (to remove various ex- 
perimental imperfections) properly. Moreover, this con- 
dition can be also profitable for an improvement of such 
experimental data for which p(p) is not reconstructed. 
For this purpose one should estimate interdependences 
between projections what is a subject of this paper. 

In the next Section we discuss the consistency condi- 
tion in its general form, i.e. for the Radon transform 
in A-dimensional space. Next, some relations between 
line (Sec. II. A) and plane projections (Sec. II. B) are de- 
rived from both the CC and a symmetry of p(p) (elec- 
tronic densisty in the momentum space) for various crys- 
tallographical structures. A fulfilment of these relations, 
proved for various model and experimental spectra, is 
discussed in Sec. III. 



II. THEORY 

The Radon transform [4] represents integrals of p(p) 
(defined in A^-dimensional space, 1Z ) over (N — 1)- 
dimensional hyperplanes: 

/oo 
p(p)S(t - C ■ P)dp, (1) 
-oo 

where £ is a unit vector in 1Z N along t and t is a per- 
pendicular distance of the hyperplane from the origin of 



the coordinate system. The equation t — £ ■ p defines the 
hyperplane. In the same coordinate system the vector p 
is described by p — \p\ and u> (e.g. Dean [5]). 

Both functions, g and p, can be expanded into spherical 
harmonics of degree I: 



s(*,C) = X>*(*)£/v(C), 



Xv 



p(p,u) = ^2pi u (p)Si u (u), 



(2) 



(3) 



Xv 



where index v distinguishes harmonics of the same order 
I. In order to make our formuale clearer, henceforth index 
v will be omitted, keeping in mind that for the same I a 
few harmonics can be used. According to Dean [5], 



Plip) = ~ 
P 



g^ +1 \t)C?(t/p)l(t/pf If-^dt, (4) 



where c = (-l) 2 ^ +1 r(; + l)r(^)/(27r^ +1 r(; + 2/i)). Here 
denotes the n th derivative of g ; Cj 1 are Gegenbauer 
polynomials and p = N/2 — 1. The equation (4) can be 
solved analitically if either [6] 



9i 



(t)^£ a;fc (l-t 2 ) A - 1/2 Q A +2fc (t), 



fc=0 



with A > N/2 - 1, or [7] 
9l(t) = e- t 



k=0 



°lkHi + 2k 



(*), 



(5) 



(0) 



where H m are the Hermite polynomials. Because in pi(p) 
all terms t™ with < n < I are equal to zero, the lowest 
term in gi is of order I. This property, called CC, follows 
from the fact that functions gi (Eq. (2)) are a linear 
combination of g{t, C,) which in turn are interdependent 
being projections of the same density. Generally, one can 
write 



9l {t)=J2 C ^ W l+2k(t), 
k=0 



(7) 
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where W n denotes an arbitrary orthogonal polynomial. 
Below, considering the cases N — 2 (line projections) 
and N = 3 (plane projections), we show that due to this 
property of <?; (t) (its minimal number of zeros) we can get 
some interdependences between projections g(t,£). For 
that it is necessary to expand them into a series of the 
same (as in Eq. (7)) polynomials: 

oo 

s(t,C) = $> fc (C)W2fc(t). (8) 

fc=0 

By combining Eqs. (2), (7) and (8) we can obtain some 
relations between c k (Q as a function of £. In the pa- 
per, for calculating the coefficients c k , the Chebyshev 
polynomials and the Gaussian quadrature formulae were 
applied. 

A. Line projections 

In the case of reconstructing p(p) from its line pro- 
jections g(x,y) (measured in the apparatus systems 
(x,y, z)), the reconstruction of the three-dimensional 
(3D) density can be reduced to a set of reconstructions 
of 2D densities (the Radon transform for N — 2) per- 
formed, independently, on succeeding planes y — const., 
parallel each other. In order to describe all projections 
in the same coordinate system (in which the symmetry 
of both p(p) and measured spectra is defined) the func- 
tion g can be characterized in the polar system where 
g(x, y) = g(t, a). Here t and a denote the distance of the 
integration line from the origin of the coordinate system 
and its angle with respect to a chosen axis, respectively. 

For the planes y = const., perpendicular to the axis of 
the crystal rotation of the order \G\, Eq. (2) reduces to 
the cosine series [8]: 

oo 

g(t,a) =J2gi(t)cos(la), (9) 
I 

with I = i ■ \G\ (1 = 0, \G\, 2\G\, etc.). Here the an- 
gle a, describing nonequivalent directions, is changing 
between < a < a G — tt/\G\. Of course, to use the 
symmetry most profitably, a sensible choice is to select 
the planes perpendicular to the main axis of the crystal 
rotation (|G| = 4 or 6 for cubic and tetragonal or hexag- 
onal structures, respectively) where a G is minimal, i.e. 
the number of equivalent directions is maximal. 

Knowing that 1/2 first coefficients cik are equal to zero 
(Eq. (7)), the equations (8) and (9) give the following 
dependences between c k ((x) 

1°. Cfc (a) = c k, i.e. \G\/2 first coefficients c k (a) are the 
same for all projections. 

2°. 2c k {a G /2) = c k {a) + c k (a G - a) for k < \G\. 

3°. 2{c k (a G /A) + c k {3a G /4)} = c k (a) + c k (a G /2 - a) + 
c k (a G /2 + a)+ c k (a G - a) for k < 2\G\. 



All these conditions have been proved for both various 
models of p(p) and for \G\ — 2, 4 and 6. 



B. Plane projections 

Due to the symmetry of p(p), in the case oi N = 3, 
g(t, £) = g(t, Q, (p) can be expanded into the lattice 
harmonics Fi(Q,tp) which form an orthogonal set of 
linear combinations of the spherical harmonics K™ = 
c(l,m)p\ m \cos<d)(e imv + e- im v). c(l,m) is the normal- 
ization constant and P ; ' m '(cos0) denotes the associated 
Legendre polynomial. Angles (<3,<p) describe the az- 
imuthal and the polar angles of the £-axis with respect 
to the reciprocal lattice system. 

In the case of the hexagonal structure Fi = K™ with 
I = mod 2 and m = mod 6 (K$ = 1, K$, K\, K%, 
K%,Kl,K%, etc.), for the tetragonal structure F t = K™ 
with I = mod 2 and m = mod 4 and for the cubic 
structures Fi are the linear combinations of K™ , where 
the first three are equal to [9]: 

F = I, 

F 4 = 0.763762611^ + 0.64549722X1, 

F 6 = 0.35355338^ - 0.9354134X|. 

So, for the hexagonal and tetragonal structures a few 
first lattice harmonics do not depend on ip. Moreover, 
very often, for some paticular sets of measured spectra, 
we cannot calculate functions gi(t) (the matrix in Eq. (2) 
is singular). For example: if for the hexagonal structure 
we have only three projections with £ along the TM, 
TK and r^4 symmetry directions (Fig. 1), the angles 
(0,(p) are equal to: (7r/2,7r/6), (tt/2,0) and (0,0), re- 
spectively. Thus, the first three lattice harmonics have 
the same values for the TM and TK directions and there 
is no a possibility to evaluate functions gi(t). Next, if e.g. 
(B,p) = (0,0), (tt/4,0), (tt/2,0), (7r/4,7r/6), (n/2,n/6), 
the set of equations for gi does not have a solution in 
spite of the fact that the first five harmonics distinguish 
here all directions £. Due to these reasons we propose 
to treat plane projections of p(p) as line projections of 
Pl and to use the consistency and symmetry conditions 
derived in Sec. II. A. 

We consider only those spectra g(t, £) for which £ is 
changed on the plane perpendicular to the axis of the 
crystal rotation of the order \G\. In such a case pL de- 
notes the line integral of p(p) along lines L parallel to 
this axis and g(t,C) can be described by Eq.(9). Now 
all relations obtained for the line projections are valid 
with c k (a) being replaced by c k (C) = c k (Q,ip) where ei- 
ther a — (p and O = 7r/2 (for the planes perpendicular 
to the main axis of the crystal rotation) or a = and 
ip = const, (for the planes with \G\ — 2, perpendicular 
to the previous ones). 

For example, for three spectra having the hexago- 
nal symmetry and with £ along the TM, TK and TA 
symmetry directions (denoted here by g(TM), g(TK) 
and g(TA)) we have the following relations: c k (TM) — 
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c k {YK) for fc = 0, 1, 2 (here \G\ = 6) and c {TA) = 
Cq(TM) = co(TK) (\G\ = 2), where Cq denotes the norm 
of the spectrum. So, to get more information about 
g(TA) we should have at least one additional spectrum 
g(TI) (for 6 = 7r/4 and e.g. ip = 0) on the plane 
with \G\ = 2 where 2ci(I7) = a(TA) + c x {YK) (see 
Fig. 1). This equality is satisfied for 6 = tt/4 and any 
if £ (0,7r/6), i.e. 2ci(r/") = Ci(rA) + ci(IW), where 
the TiV and TJ" directions are defined by the same ip 
with 6 = 7r/2 and 6 = tt/4, respectively. Knowing that 
Ci(IW) = ci (does not depend on (ys) we obtain that 
ci(0 = 7r/4, 93) has the same value for each spectrum 
g(C) with £ described by = tt/A and any <p. This is 
derived by treating g(TN) simultaneously as the line pro- 
jection of p|6| and p|2 with the symmetry |G| = 6 and 
\G\ = 2, respectively. Here p| 6 | represents the line projec- 
tion of p(p) along the main axis of the rotation, while p~\2\ 
along any line perpendicular to this axis. This last depen- 
dence is particularly interesting because it gives interde- 
pendence between line projections of different densities, 
while all consistency conditions for the line projections 
are for the same p. 

All dependencies, shown here on the example of the 
hexagonal structure, are valid for the cubic and tetrag- 
onal structures where \G\ = 6 is replaced by \G\ = 4. 
However, for the cubic structures, where three axes of 
the fourth order exist, one can get some additional rules. 
Because directions (6, ip) = (0, 0) and (7r/2, 0) are equiv- 
alent, we obtain that not only c but also c\ must be the 
same for all projections. Some of these results can be also 
derived from Eq. (2) because for the cubic structures all 
lattice harmonics (except for F ) depend on (0, p). Tak- 
ing as an example directions [OOh], [klO] and [klh] we 
obtain the following dependences between a: 

oiCi([00ft]) + a 2 c t ([kl0]) + a z c,{[klh]) = for i = 0, 1 
hcidOOh]) + b 2 Ci([klO]) + b 3 c t ([klh}) = for i = 0, 1, 2 

where J2 &» = J2 a i = °- 

The above equalities arise from the conditions of van- 
ishing all polynomials Wi with i = 0,1 and i = 0,1,2 
in the expansion of g± (t) and ge(t), respectively. The 
coefficient c\ satisfies both equations when ci([00ft.]) = 
Ci([fcZ0]) = c\{[klh\) what is in agreement with the previ- 
ous result. Knowing additionally that Ci([fcZ0]) does not 
depend on the direction [klO], we can conclude that for 
the cubic structure the value of c\([klh}) does not depend 
on the direction (denoted here by [klh]). 

As before, all relations were proved for both model 
and experimental profiles. Some examples are presented 
in the next section. 



III. APPLICATION 

First, disposing twenty five 2D ACAR experimental 
spectra [10] for the four metals of the hep structure with 
£ on the plane \G\ = 6, we obtained that the conditions 1° 
and 2° are satisfied with a very high accuracy, lower than 
0.5% (ci), 1.5% (c 2 ), 1.5% (cs), 2.5% (c 4 ) and 5.5% (c 5 ). 



Because the first ckm(Q) have the highest values, they are 
determined the best (the influence of the statistical noise 
is the lowest). This very small inconsistency implies that 
spectra [10] (with the average number of counts at peak 
about 60000) were not only measured but also corrected 
with a very high precision. 

In the figure 2 we show results for the theoretical 
Compton profiles for Cr [11] with £ along [100], [110] 
and [111] directions. In this case (cubic structure) the 
first two coefficients Ck{C) should have the same value 
for each spectrum, while the third coefficient should sat- 
isfy the relation: c[100] - c[110] = 3{c[110] - c[lll]}. 
After normalizing spectra to the same area (co = 1) we 
obtained that a distortion of c\ from its average value 
is changed from 0.1% up to 0.4% and the inconsistency 
of C2 is of order 1.5%. Thus, the behaviour of Ck(C) is 
the same as for previously studied experimental data. It 
is connected with the fact that in order to create such 
spectra it would be necessary to calculate p(p) in the 
whole momentum space p, what, of course, cannot be 
done. Here the accuracy of calculating Compton profiles 
was comparable with the statistical noise of 2D ACAR 
data [10]. Of course, in order to get the similar inconsis- 
tency of the Compton profiles (where a contribution of 
core densities is much higher than in a positron annihila- 
tion experiment) the experimental statistic must be also 
much higher. 

In order to examine how these conditions react to 
an improper shape of g(t,C), we changed somewhat the 
shape of one spectrum ( marked in figure 2 by squares). 
Here we would like to point out that this incorrect spec- 
trum does not differ from other spectra too much (the 
differences between g([100]) and g([lll]) are higher). Af- 
ter normalizing it to the same area (co = 1) we obtained 
that Ci changed its value from -0.7 to -0.74 . So, here 
we observe a much higher inconsistency (5%) than for 
experimental data [10] (0.5%). 



IV. CONCLUSIONS 

All spectra received from experiments are contami- 
nated by statistical noise and thus they are not consis- 
tent. Depending on the number of measured projections, 
some part of the inconsistent noise is eliminated via the 
consistency conditions during the reconstruction of p(p). 
However, we propose to check (before the reconstruction) 
if spectra were measured and corrected properly (the in- 
consistent part of the data cannot be too large). For 
that purpose we can use the interdependences, derived 
by combining the consistency and symmetry conditions, 
obtained here for the line and plane projections. 

The knowledge of the interdependences rules can be 
also profitable for an improvement of data for which den- 
sities are not reconstructed (some part of the inconsistent 
noise can be eliminated). Moreover, it can be applied 
to verify these techniques which are used for correcting 
Compton scattering spectra (they are not univocal and 
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can be individual for each spectrum [12]). For that we 
propose to measure the high resolution CP, denoted usu- 
ally by J(pz), for hexagonal metals with p z = £ changed 
on the plane \G\ =6 (6 = 7r/2). Having p z along TK 
(if = 0), TM {ip = tt/6) and TN (ip = tt/12), we can 
check if a(TK) = a(TM) = a(TN) for i = 0,1,2 and 
if Ci(TK) + Ci{TM) = 2ci(IW) for i = 3,4, 5. Of course, 
the best choice is to study rare-earth metals where the 
anisotropy of p(p) is so high that spectra J(p z ) should 
be essentionally different. 
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FIG. 1: Planes and directions for the hexagonal structure 
having symmetry of the sixth and the second order, marked 
by dots and lines, respectively. 




FIG. 2: Theoretical plane projections for Cr [11] having cu- 
bic symmetry with £ along [100], [110] and [111] directions 
marked by full, broken and dotted lines, respectivelly. Incor- 
rect projection [110] is marked by the full squares. 



